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T— I Abstract 

o 

fvq In this paper, we propose a domain decomposition method for niuhiscale second order eUiptic 

» partial differential equations with highly varying coefficients. The method is based on a discon- 

^ C^ tinuous Galerkin formulation. We present both a nonoverlapping and an overlapping version of 

^> the method. We prove that the condition number bound of the preconditioned algebraic system 

in either case can be made independent of the coefficients under certain assumptions. Also, in 
our analysis, we do not need to assume that the coefficients are continuous across the coarse 
grid boundaries. The analysis and the condition number bounds are new, and contribute to- 

I— I wards further extension of the theory for the discontinuous Galerkin discretization for multiscale 

•^ problems. 

^ 1 Introduction 

a 

2 Subsurface flows in heterogenous media [Sllin], that is, where the heterogeneity varies over a wide 

'""' range of scales, are examples of multiscale problems. Numerical solutions of such problems are often 

,—1 affected by the heterogeneity, in particular when it is highly varying, that is, the permeability of 

K*" the media can span a large scale. In this paper we consider the numerical solution of flow problems 

^L governed by elliptic equations with highly varying coefficients. Domain decomposition methods use 

(^ both fine scale and coarse scale subproblems as well as interpolation operators from the subspaces to 

'^ the solution space to construct preconditioners for the original problem. The convergence property 

cn is linked to proper treatment of the jumps in the coefficients [5]. A key point when building such 

Z^ domain decomposition methods is to find a good coarse problem that can capture relevant small 

^— I scale information from the fine level. In general, due to the complex geometry of the conductivity 

^ field, where high and low conductivity regions often appear as small inclusions inside subdomains 

and long channels across subdomain boundaries, it is rather difficult to design a robust domain 
decomposition method for such problems. 
C^ When considering high conductivity regions (inclusions or channels) in a low conductivity back- 

ground, the authors of [5j used functions that are discrete harmonic in each coarse grid block with 
special boundary values as coarse basis functions (multiscale basis functions). They introduced 
two indicators, 7r(a) and 7(a), where a is the coefficient representing the conductivity, to capture 
the effect of the jumps in a on the condition number. The indicator 7r(a) measures how well the 
coarse partitioning of the whole domain i} is, as for instance, 7r(a) behaves well as long as the high 
conductivity regions (inclusions) do not cross any subdomain boundaries, and badly otherwise. 
The indicator 7(a) is a measure for the weighted energy norm of the coarse basis functions, and 
it depends on the choice of boundary conditions, as for instance, an oscillatory boundary condi- 
tion is often needed in order to keep the energy norm low, consequently for 7(a) to behave well. 
The authors in [5] used the traditional L^ projection operator from the fine space to the coarse 
space, and proved that the condition number bound is independent of the jumps. This approach 



X 



assumes that the conductivity coefficient is continuous across coarse grid boundaries, although 
their numerical results did not seem to require this. Recently, in [11], the authors proposed an 
overlapping domain decomposition method for the multiscale problem. The idea of their method 
is based on the fact that high conductivity regions correspond to the smallest eigenvalues of the 
system. Consequently, they used the corresponding eigenfunctions as coarse basis functions, and 
proved a weighted Poincare inequality resulting in a condition number bound independent of the 
jumps. They also proposed an overlapping domain decomposition methods for the Schur comple- 
ment system. 

Discontinuous Galerkin methods may offer several important and valuable computational ad- 
vantages over their conforming Galerkin counterparts. The finite element spaces are not subject 
to inter-element continuity conditions, and the local element spaces can be defined independently. 
This makes discontinuous Galerkin methods well suited for their applications to multiscale prob- 
lems with piecewise constant coefficients relative to a fine triangulation. A domain decomposition 
method for the discontinuous Galerkin formulation of a multiscale elliptic problem has recently 
been proposed in [1]. There a composite discontinuous Galerkin formulation, that is, a regular 
continuous formulation inside each subdomain and a discontinuous Galerkin formulation across the 
subdomain boundaries, has been used. There the coarse space consists of piecewise constant basis 
functions over the coarse partition. The condition number bound of this method is shown to be 

ai Hi 
maxmax — — — — , (1-lj 

i j Oj hihij 

where Oj and a^ represent the maximum and the minimum of the coefficients inside a boundary layer 
of the subdomain fij. Hi is the diameter of Oj, hi is the fine mesh size in Oj, and hij = 2hihj/{hi+hj), 
the harmonic average of hi and hj. We note that the condition number bound above depends on 
the jump of the coefficients inside the boundary layer of each subdomain. 

The present work is an extension of the work in [5] to a discontinuous Galerkin formulation. 
We use the same bilinear form as the one given in [13] with harmonic average weight functions, 
defined on the fine space which is the space of piecewise linear polynomials with respect to the 
fine triangulation. A composite discontinuous Galerkin formulation, c.f., |4J, is used on the coarse 
space whose basis functions are the multiscale basis functions given by the oscillatory boundary 
conditions, c.f., [5]. We present both a nonoverlapping and an overlapping method. 

A new indicator /3(a) is introduced measuring the L^ norm of the jump of the coarse basis 
functions across coarse grid boundaries (see (3.11|)). We show that, under certain assumptions. 



slightly weaker than in [5|, our methods are robust. For the nonoverlapping case, we show that the 
condition number bound is 

77 max max We^{l)-, — h max (7(a), r//3(a)), (1-2) 



where We is the harmonic average weight function (see (2.3)), 7(a) and /3{a) are the indicators, 
and rj is a penalty parameter. The numerical experiments support the assumption that this bound 
is sharp. However, due to the presence of We in the bound, the method may have a large condition 
number when high conductivity channels cross subdomain boundaries. Note that, this is the case for 



( 1.1 ) even when there are no channels crossing subdomain boundaries. The situation gets improved 



in the overlapping case. We show that, for the overlapping method, the condition number bound 



IS 

TJ , 

■K^a) max (7(1), /3(1)) max — ^ + max (7(a),r//3(a)). (1-3) 

i Oi 

The rest of the paper is organized as the following. In Section 2, we introduce our model 
problem and the discontinuous Galerkin discrete formulation. The two level additive Schwarz 
domain decomposition methods, first the nonoverlapping and then the overlapping version, are 
defined and analyzed in Section 3. Section 4 is devoted to the numerical results. 

2 Problem Setting 

Throughout this paper, we adopt standard notations from the Lebesgue and Sobolev space theory 
(see [1]). We further use ^ < i? to denote A < CB with a positive constant C depending only on 
the shape regularity of the meshes, and A^ B to denote A< B < A. 

Consider the self adjoint elliptic problem on a polygonal domain il with boundary 9rj, 

L{u) = -V • (a(x)Vu) = /, X G f^, (2.1) 

u = g, xedn, (2.2) 

where the coefficient a{x) G L°°{Q,) with a{x) > ao > 0, representing the conductivity field on Q, is 
piecewise constant on the fine mesh Th, a{x) may have (discontinuous) jumps across the elements. 
Also we let /G L'^{n). 

First, we introduce some notations. Th denotes a fine triangulation of the whole domain O 
which is quasi-uniform, Th = UrerJ ''"' where r represents a small triangle in Th and h is the mesh 
size. Th denotes a coarse triangulation which we get by partitioning 0, into triangular substructures 
Th, which is also quasi-uniform with the mesh size H. We assume that the boundary edge of each 
coarse grid element in Th , is aligned with the edges of the elements in the fine triangulation Th ■ We 
denote the elements of the coarse and the fine triangulation by K, and r, respectively, e denotes 
an edge of a fine element r, and £ is the union of all edges in $1, <S = IJecOr ^- Additionally, 
£ = £j \J£d, where £1 := Uj>i ^ji refers to all the inner edges and £d '■= \Jecdn ^ refers to all the 
edges touching dQ. Given a coarse triangulation Th = {-f^ili^i, we let dKi be the boundary of the 
element Ki, and dKij = dKi n dKj be the open edge shared by the elements Ki and Kj. 

Next, we introduce two weight functions related to each fine edge e £ £1. We first denote the 
two fine elements sharing an edge e by r^ and r^ , and denote the coefficients of the two elements by 
a^ and a'L respectively. The weight functions w^ and w^, associated with the edge e, are defined 
as 

wl+wt = 1, 

where w^ = a!i/(a^ -|- a'L) and wL = a+/(a+ + al). The harmonic average weight function, 
associated with the edge e, is then 

20-^ o-^ 
We = wla\ + wtaL = - + . (2.3) 

a_ + ar^ 

The following inequalities hold. 



We wLaL + tfia": 



and 



W, 



< 1. 



Now, let amin = min(a!i,a^) and amax = max(a^,a^), it then follows immediately that 
Oimin ^ We < 2amin- In the above definitions of We and w'^, for e G iS/), we set ti^e = 1 and 
We = a". 

Let the jump across an edge e be 

[.] = /"+"+ + "-"-' ^^^^' (2.5) 

where n^ and n^ denote the unit outward normal vectors of r^ and rt , respectively. The weighted 
average {v}^ is defined similarly by 

M^=(->^ + --"-• ^^f- (2.6) 

For the ease of our presentation, we define all our norms and function spaces here. 
Given a domain D <Z 0,, we define the standard L°^ norm on D as 

ll«llo,oo,D := ess sup{|u(x)| : x e D}, 
and the standard L^ norm on D as 

ll^llo D •= / ^ dx. 

Jd 
A weighted H^ seminorm on D is defined by 

' ' Jd 
where a{x) > 0, Vx G D. With a{x) = 1, Vx G D, we get the standard H^ seminorm 

? n := / Vn • Vudx. 



\u\i,D ■ 

Id 



A weighted norm, based on the discontinuous Galerkin formulation, on D, is defined by 

II 1 1 2 \~~^ It—! t—! 7 \~~^ ^e 1 1 r 1 1 1 2 \~~^ " e 1 1 1 1 2 

ll^lli,h,a,D := 2^ / aVu-Vudx+ ^ -7-||Nllo,e+ 2^ 



rCD-''^ ^jlCD "= eGdD ^ 



where U7>/ ^ji refers to all the fine element edges in the interior of D. We use 

\lh,c.,D-= E favu-vudx+ Y^ ^mwi 



to denote the corresponding semi norm on D. If a{x) = 1, Vx £ D, we get the fohowing norm 



\lh,D-=Yl ['^U-Vudx+ Y^ ^^\[u]\\le+ Yl /T 






\lh,D-=Yl /Vu-Vndx+ Y, jr 



and seminorm 

\u\lh,D-=yZ I Vu-Vudx+ Y 7-|IMIIo,e- 

j>l 
Let H^{Th) be a broken Sobolev space of degree s > defined as 

H'iTh) = {ve L\n) : v\r G H'{t), Vr G %}, 

and F*(7h) be its subspace defined as 

V%Th) = {vG H^Th) : V • [Vv) e L\t), Vt G TJ. 

With the above preparation, we can now define our discontinuous Galerkin bihnear form a{u, v) 



and the right hand side f{v) for the continuous problem (2.1) and (2.2). They are defined as 
follows, c.f., [13J. For u,v £ V^^'{Th), 



i{u,v) =2^ / aVu-Vvdx — 2, I {(ySJu}l^[v\ds — 2, / {aV v}l^[u]ds 

+ V Y IT MM'^*- Y WeOeVUe ■ Velleds (2.7) 

Y / W^eOeVVe • UeTleds + V Y JT '^eVeds, 



ecan"^ ecan 



and 



where Uj^iCji refers to all the fine element edges in the interior of O, Wg is the harmonic weight 



function defined in (2.3), and r/ is a penalty parameter. 



Here and in the text below, for simplicity, we use We, he and {■}% instead of VFe, , he i and {-^w 



respectively if there is no confusion. Moreover, as stated in |13j . for / G -L^(O), the flux a = —aVu 
is in H{div;n) D H^Thf , where iJ(div; fi) = {v G L'^{^f : V • v G L'^{9)] is an Hilbert space 
equipped with the norm ||v||//^.^.j5 = (||v||q ^ + ||V • v||q j^)^/^. The integrations on the fine edge e^/ 



in (2.7) above can be understood in the weak sense, see [13] for details. 

Due to [13], the the weak solution of (2.1) and (2.2), u G H^+'{^) n V^+'{Th), < e < s < 1, 
satisfies the following variational equation 

a{u,v) = f{v),yveV^+\Th). (2.8) 

Next, we define the finite element space V^ associated with Th as follows. For any fine triangle 
r let -Pi(t) denote the set of all linear polynomials on r. The finite element space V is then defined 



as 

V^ = {v:v\^ePi{T),yT(^Th}. 

The bilinear form on the finite element space V is defined as ah{u,v) = a{u,v), \/u,v € V . 
We can now formulate our discrete problem: Find Uh S V^ such that 

ah{uh,Vh) = f{vh),yvh€V\ (2.9) 

Naturally, the above bilinear form induces a norm in the space V^, which is 

II ii2 \~^ It—? t—?i \~^ ^Ciir 111 2 \~^ '^e II ii2 

\Ml,h,a,n = Z^ aVv-Vvdx+ ^ ^IIMIIo,e + 2^ -^|Pe|lo,e- 

3>l 

The next lemma gives us the continuity and coercivity of the bilinear form a/i(-, •). 
Lemma 2.1. There exists a constant C, such that 

(i) \<^h{u,v)\ < cMli^^^^^Mlf^^^^^, yu,v G VK 

(ii) There exists positive constant r]Q > 0, such that for all rj > rjQ, 
where C{rf) is independent of the jump in the coefficient. 



As a consequence of the above lemma, (2.9) has a unique solution. The proof of Lemma 2.1 
well as an error estimate is given in [13j . 



as 



3 The Schwarz methods 

In this section, we use the Schwarz framework [21 [12] to design and analyze our additive Schwarz 
domain decomposition methods for the discontinuous Galerkin formulation. Let J7 be partitioned 
into a family of nonoverlapping open subdomains {r^j, 1 < i < N} with 

N 

= [J ni, rjj n Oj = 0, i 7^ j. 

Next, we give the overlapping partition {Q.'j}f^^ by extending each subregion 0.1 into a larger region 
rj'j, i.e., Oi C O'j, so that 

N 

n = U o[. 

i=\ 

We assume that both the nonoverlapping and overlapping partition are aligned with Th and Th- 
We denote by 5 the minimum of the distance between the boundaries of Oj and il^, i.e., b = 
min dist{dQ'^\dO,dQi\dO,). If there exists a constant number C > such that S > CH, we say 
that {O^}^^ has generous overlap, and if 6 is proportional to h, we say it has small overlap. 

Following the Schwarz framework, the space V^ is split into a number of local subspaces and 
a global coarse space, i.e., V^ = Vb + J2i=i^i^^ where Vq is the coarse space and {V/^}^i are 
the local subspaces. For the coarse space, Vq, we use the standard multiscale finite element basis 



functions as described in [SjIS], where, we allow the basis functions to be continuous inside each 
coarse element K, and discontinuous across dK. Below is the description of Vq. 
We define our coarse space associated with Th as follows 



Vo = span{(l)p,K, for all Xp G M^ (K), VK G Th}, (3.1) 



where M (K) includes all the vertex nodes of K, and (jyp^K denotes the multiscale basis function 



defined in (3.3). To be more specific, we need to introduce suitable boundary data ipp^oK, which is 
required to be piecewise linear(w.r.t. the given fine mesh 7h restricted to dK), and to satisfy the 
following: %,dK{x'j^) = ^p,p', V Xp,x'p G M^{K), < ipp,dK{x) < 1, andJ2xp(^AfH(K) ^P,dK{x) = 
1, V x G dK, and ipp^dxix) = 0, on the edge of K opposite Xp. 

An obvious choice for the boundary data i^p^oK satisfying the above conditions is the standard 
linear boundary condition. 

The linear boundary condition works well when the high conductivity regions lie strictly inside 
coarse grid blocks. However, when they touch the the coarse grid boundaries, linear boundary 
condition fails. In this case, we use another boundary condition, also known as oscillatory boundary 
condition, c.f., f5], which satisfies the above conditions and is effective. The description of this 
boundary condition is as follows: 

Let T be an edge of the coarse mesh Th with end points Xp and Xp/ , and a be the restriction 
of a to T. Then the oscillatory boundary condition is given by the finite element solution of the 
following two-point boundary value problem: 

-{a'^iijjyy = 0, X G T 

ipjixp) = 1, and ipj{xp>) = 0. 

Since the coefficient a is piecewise constant, the finite element solution of the above equation 
can be expressed explicitly by 

^J(x) = ( / {a^)-^ds)-\ I {a^)-^ds), for ah x G T, (3.2) 

where T^; denotes the line from Xpi to x. The function ■0J is continuous and piecewise linear with 
respect to T CiTh- We set ijJp^dK = V'J on each edge T of K, containing Xp, and V'p,(9-ft:|T = on the 
edge opposite to Xp. 

Once the boundary condition tpp^oK is determined, (pp^K is constructed by a discrete harmonic 
extension inside K. First, define the Pi-conforming finite element space associated with Th as 

S (Q) = {v : v\r is linear and u is a continous function over Q,}. 

Then (f)p^K G S^{K) can be defined as 

/ aV4>p,K ■ Vvh = 0, for all vt G So{K) subject to (l)p,K\dK = ipp,dK, (3.3) 

JK 

where Sq{K) = S^{K) n Hq{K). The most important property of (/)p^K is the energy minimizing, 
which can be stated as follows, 

\4>p,K\i,a,K < \()\i,a,K, for ah G S'^iK) which satisfies 8\qk = 4'p,K\dK- (3.4) 



Note that our coarse space Vq actually includes functions which are required to be continuous 
inside each coarse element K and discontinuous across coarse grid boundaries, and obviously we 
have Fo C y'^. 

Having the above preparations, we may define ao(-, •) the bilinear form associated with Vq as 



ao{u,v) =y^ / aVu -Vvclx — N^ N^ / {aVM}^[f](is 

rn>k j^l 

m>k j^i m>k j^i yO.O j 

- ^ ^ / WeOeVUe ■ Velleds - ^ ^ WeU^Ve ■ Uelleds 
+ ^ E E J-^'^eVeds, yu,V eVo. 

Note that Vq C F^, and ao(u, v) = ah{u,v), \/u,v G Vq, hence it follows from Lemma 2.1 that 
ao(v,v) is coercive on Vq- 

Next, we define the local spaces, {Vl^}f^i associated with subdomain partition {i?j}^i. 

Vt = {veV^ ■.v = {)\n n\Bl}, i = l...N, (3.6) 

where Bi = fij for the nonoverlapping partition, and Bi = Q'- for the overlapping partition. 
The corresponding local bilinear forms ai{u,v) can be defined as follows: 

ai(u, f) = Y^ aVu-Vvdx- V^ / {aVn}^[u](is — Y^ / {aVv}l^[u\ds 

3>l 3>l 

~^^ E / 1~MN^S- E / WetteVUe ■ VeUeds 

— 2, / W^eOeVlie • Ve^Veds — N^ / WeOLeVVe ■ UeT^eds 

+ "?/ 



We. 



-UgVeds, Vu, f G VJ , 



/e ""£ 



2.1 



that 



where T^ = dBi\dn, H^ = 95^ n dn. 

Note that V^'^ C F'*, and ai{u,v) = ah{u,v), \/u,v G V^'^, hence it follows from Lemma 
ai{v,v) is coercive on Vj^. 

Remark 3.1. Note that the local bilinear form aj(-,-) on V/^ x V^^ is actually the restriction of 
a/i(-, •) on Vf' X Vf^. Thus, in the implementation of this method, after we build the global stiffness 
matrix A corresponding to the bilinear form ah{-,-), we can easily get the local stiffness matrix Ai 
corresponding to ai(-, •) by taking the i-th diagonal block of A. 

With this preparation, the two level additive Schwarz domain decomposition method can be 



presented as follows. 

For 1 < i < A^, we define the operators Tj : V^ — )• V/^ by 

ai{TiU,v) = ahiu,v), Vn gV'', v£ VJ", 

and for i = 0, we define the operator To by 

ao(Tou,v) = ahiu,v), Vn G V^, v G Vq. 

Clearly, each of these problems has a unique solution. 
We define the additive operator as 

T = To + Ti + • • • + Ttv, 



and replace (2.9) by the operator equation 

Tuh = fh, (3.8) 

where 

N 
fh = '^gi with gi = TiUh, 



i=0 



and Uh is the solution of (2.9). 



3.1 Analysis 

In this section, we estimate the condition number of both the nonoverlapping and overlapping 
additive Schwarz method. We use the standard Schwarz framework [2lll2j. 

We will see that, for the nonoverlapping method, the condition number will depend on the two 
indicators 7(a) and /3(a), while for the overlapping method, it depends on an additional indicator 
7r(a) which will be defined in Section 3.1.2. 

The first indicator 7(a), which is borrowed from [5], measures the maximum weighted energy of 
all the coarse basis functions, which can be used to have an indication of how well the coarse basis 
functions are constructed. We will show that, by choosing suitable boundary conditions, 7(a) can 
be bounded independently of the jumps in the coefficients. 

Coarse robustness indicator [5|. Given a coarse triangulation Tj/, and the set of coarse basis 
functions {4>p^K, Vxp G M^{K), \/K G Tff,} then 

7(a) = ma^ max \(t)p,K\l,a,K- 

Next, we introduce the indicator /3(a), which measures the weighted L2 norm of the jump of the 
multiscale basis functions on the coarse grid boundaries. It is the maximum value of /3 (a), which 
corresponds to the integration on the inner coarse grid boundaries that are inside il, and /3^(a), 
which corresponds to the integration on the coarse grid boundaries which intersect with dO,. The 
term /3(a) enters into our analysis due to the use of the DG bilinear form on the coarse space. 

Multiscale DG indicator. Given a coarse triangulation Th, and the set of coarse basis func- 



tions {(j)p^K,^Xp G M^{K), ^K eTn,} then 



(3\a)= max max V ( / ^UpM'^ds). (3.9) 

where [0p] = 4'p,Km ~ 4>p,Kk ■• dK^k is the edge shared by the coarse elements Km and Kk , We is 



defined in (2.3), and h^ is the length of the edge Cji. 



However, when dKm H dil / 0, we set 

1^^^^^ = ^^^\c. ™f^ X] i ir\(l^p,mfds). (3.10) 






Having all the above preparations, we define 

/3(Q) = max(/3-^(a),^-^(a)). (3.11) 

3.1.1 Nonover lapping additive Schwarz method 

In this section, we propose a two level additive Schwarz method with nonoverlapping subdomains, 
and present an analysis of the condition number. Since it is a two level additive Schwarz method, 
we use Vq in the previous section as our coarse space with bilinear form ao(-,-)- ^'^^ the local 
subspaces, by taking Bi = Oj in (|3.6l), we have 



vj" = {vev'' :v = om n\ni}, i = 1 . . . iV. 



Note that, in this case, V is a direct sum of {V^ }• The local forms aj(-,-), 1 < i < A^ are 



given by (3.7) with Bi = Jlj. 

In order to estimate the condition number of the two level nonoverlapping domain decomposition 
method, we need to define an interface bilinear form, /(•, •) : V x V — )• R, as follows: Vu, v G V , 



I{u,v) = - 

aK. 



y ^ y { / ajWjVuj ■ vinids + / ajWjVvj ■ uinids 

m>k jyl 

+ / aiwiVui ■ Vjiijds + / aiwiVvi ■ Ujiijds \'^-^^) 

+1] -j-^UjVids + 1] —^uiVjds}, 

Jeji "e Jeji "e 

where Uj = u\r- and ui = ulyp with Cji being the common edge between tj and r/. 

Observe that the relationship between the bilinear form on the fine space and the bilinear forms 
on the local subspaces is given by 

N 

ah{u, v) = ^ ai{ui, Vi) + I{u, v) Vu, v G V^ , m, Vi G Vj^, (3.13) 

where u = X]j=i ■"«' ^^^ ^ = ^j=i Vi. 



For theproof of (3.13), one only need to compare the terms in ah{u,v) with those in J2i=i (^i{ui,Vi), 



10 



The next lemma which states the Poincare and a trace inequahty for our discontinuous case, 
has been proved in J14j . 



Lemma 3.1. Let D he a convex domain, {Tb} he a family of partitions of D, and D = UtsT ^ 
with diam{T) ss h, then for any u ^ Vd = lircTo^^i'^)-: ^^^^ ^ = \h\ Id"^^^ being the average 
value of u over D, we have 

\W - u\\o,D < diam{D)\u\i^h,D, 



where 



\lh,D=Y.\\^<r+ E ^"'IIMIIL 



tCD EjiC-C 



and [u] denote the jump on the edge eji. Consequently, if diam{D) ~ 0{H), we have 

\MldD ^ H^^IMId + H\u\lu n. 



Remark 3.2. Note that, in the above Lemma, the condition on the convexity of the domain D is 
actually too strict. It is stated in p^ that this can be dropped. 

Define the restriction operator Rh '■ V (Q,) — )■ Vo{Q,) as follows: 

Rhu{xp) = Tip, Vn G V^\n), Vxp G N"{K), VK G Th, 

where Up = rj-r f^ udx, with u!p being the union of elements sharing Xp, as shown in the figure 
below. 




Lemma 3.2. For the restriction operator defined ahove, the following approximation and stability 
properties hold 



\u - RhuWIk < CH'\u\lh,^^^, Vn G V^Q). 



ao{RHU,RHu) < C max [j (a), r]j3 (a)) a h{u,u), Vu G F (il). 
Particularly, if a{x) = 1, \/x £ fi, and rj = 1, then we have 

aoiRHU,RHu) < Cmax(7(l),/3(l))a/j(n,n), Vu G V^{n), 

where ujk = \Jp(^j\fH(x) ^pi ^^^ ^^^ constant C is independent of h and H. 
Proof. Since 

Up = -. r / udx < C- r||-u||o,ajp|Wp|2 < C rikllo.wp, (3-14) 



11 



we have 



\u - Rhu\\Ik <C||ii|lo,x + C\\Rhu\\1^k 






\uf 



XpGAfHiK) 



where in the last inequahty we have used \ujp\ > \K\ and Up is part of lok- 

Since our interpolation operator Rh keeps constants unchanged, we can take u — j^^ f^ udx 

instead of u, and the approximation property holds due to the Poincare inequality in Lemma |3.1[ 
Next, we prove the stability property. Let 

ao{RHU, Rhu) = y^ / oSJRhu ■ VRnudx 
K -^ 

-2 Yl Y. I {(^'^RHu]l[RHu]ds 

m>k jyl 

+ '? E T. [ ^[RHu][RHu]ds 

m>k jyl 

- 2 ^ ^ / WeaeV{RHu)e ■ {RHu)e^eds 
+ ?? ^ ^ j-^{RHu)e{RHu)eds 



dKcOn eCdK 



=h-2l2 + h-2h + h. 
We begin by estimating the term Ii. For each K G 7h, with u = u — -rj-^ f^ udx, it follows 



from the Poincare inequality in Lemma 3.1 that 



\Rhu\I^^K = \RHu\la,K ^ Yl (l^pl ^PIIo,a;p)l'^P.^Il,a,i^ 

XpeM"{K) 

Summing over the elements K G Th, we get 



h ^ 7(a) Yl \^\lh,i^K - l{a)ah{u, u)- (3.15) 



K 
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Next, we estimate the term J2. Let v = Rhu, then by the Cauchy-Schwarz inequahty 






^ E E (/ ^|{"V.};li^^«)+ E E / f^iHi'^- 



m>k jyl 






Prom the definition of {-j^, and (2.4), we have 

i^2i< E E (/ ^K«iv^i+<«;v<|2d.) 

aK_i.cn e„-,caK_i. ■'Sji "^ 



m>k jyl 

+ E E 



■IHPd. 



m>k jyl 



< 



E E 



/le max( 



(«;^)2al («^^)2a^,, „,^^„,2 , „e,v7..e,2^ 






+ E E 



W, 



w. 



'-){a'L\Vv'L\' + al\Vvl\')ds) 



W, 



MPds 






< 



E E 






+ E E 



hea'L\Vvl\'^ds+ I hea'+\Vvl\'^ds) 



en 



TV. 



M?ds. 






Let T^ and r.f be the two fine elements sharing the edge e, and note that Vv!j_ and Vvi. are 
constants on t1_ and r.^ respectively. Then 



^2|< E E {Sa-\^V-?ds+fa\\Vvl\''ds) 



m>fe j>; 

+ E E 



TV, 



IHpd. 






-M?ds 



^ E E (/ j«^/^v.pdx)+ Y. 1: I 

SK^kCn e^iCdK^k JT^VJrt SA'^feCn ejiCSA'^fc "^e^ 



(3.16) 
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We now consider the term I3. First we note that, for ah dKmk, m > k, we have 






e he 



\[RHu]\'ds 



CdK 



jl'^S'^mk XpeN'H {Km) '' ^^^ ^ 






U„ 



Xp€J^H(Km) 



oeaK„ 



oedK^ 






\0,^K„ 



XpeMH[Km) 



ieaK„ 



Taking u = u— -. — ^ — , f udx instead of u in the above equation, it foUows from the Poincare 
inequahty in Lemma |3.1| that 

Y j ^\[RHufds=^ Y: I ^mHu]?ds 

<r?/3^(a)|n|2. 






<^jl<^SKmk ^^J' 






,^Kr, 



Hence, we have 



^3 = .E E 



We 
he 



9Kmk e.iCdKmk '''' "" 

?n>fc 



[RHu]\'^ds < ^/3^(")El^ll^.'^if™ ~ vP'{<^)ah{u,u). 



(3.17) 



Finally, for the terms on the boundary dQ, we apply the same techniques, and we get 

|/4|</l+/5, (3.18) 



and 



^5 = ^ E E /^(^^")-(^^")-'^"^^/5''(")Ei^i?A 



^Kr, 



dKcdfleCdK 



(3.19) 



<7 



ri(3 {a)ah{u,u). 



The stability estimate for the general case thus follows from (3.15), ( 3.16| ), (3.17), (3.18) and 



(3.19) above. The estimate for the particular case holds naturally. 



We can now give an explicit bound for the condition number of our two level nonoverlapping 
additive Schwarz method. 

Theorem 3.1. For all u G V^ , there exists uq G Vq and Ui G V/^, I < i < N, such that 

N 

Yo'ii^i^'^i) ^ CXah{u,u). 

1=0 

Consequently, we have k.{T) < A, where k{T) denotes the condition number of the additive Schwarz 
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operator T as defined in (3.8), and A is given as 



A = r/ max max VFe7(l)— + max (7(Q),ry/3(a)), 

m>fc 

where C is a constant independent of h, H, and a. 

Proof. We need to verify the three assumptions of the Schwarz framework |2l [12] . More precisely, we 
need to estimate the three constants Cq, p{M) and w which corresponds to the three assumptions. 
The bound of the condition number is then given as: 

K{T)<Cluj{p{M) + l). 

The first assumption of the Schwarz framework, asks for an estimate of the smahest Cq such 
that 

N 



^ ai{ui, Ui) < Clah{u, u). 



i=0 
rh 



For all u G V , let uq = Rhu, z = u — uq and Ui = z\n^. From (3.13), we know that 

N 
ah{u -Uo,U- Uq) = ^ ai{ui, Ui) + I{u -Uq,U- Uq). 

4 = 1 

Hence, 

N 

\'Y^ai{ui,Ui)\ < \ah{u-uo,u-uo) - I{u - uo,u - uo)\ 



i=l 



< \ah{u-uo,u-uo)\ + \I{u-uo,u-uo)\ = h + h- 



m>k j^l 



(3.20) 



Using the fact that ah{uQ,UQ) = ao{uQ,uo) and Lemma 3.2 it then follows 

h = \ahiu-uo,u-uo)\ < ah{u,u) +ao{uo,uo) < max (j (a), r] (3(a)) ah{u,u). (3.21) 

Next, we estimate the term 12- 

h = \I(.z,z)\ < ^ ^ I / "J'^J^^J' ■ ^'"'"^^l 

m>k j^l 

m>k jyl 

+ ^ X] X^ I / -ff-Zj^ids\ = hi + h2 + -^23- 
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To estimate the terms I21 and I22, we use similar techniques as in Lemma 3.2, which gives 

7^ air ^ ,. ,-flf _ "-e Jen 



'212, 






and 



h2<Y.\\a'/'Vz\\l^+Y: E ^if \^j\'ds) = h2i + l222. 



m>k j-^l 



Using the Cauchy-Schwarz inequahty for the term I23, we have 

We 



123 =v Y. E 



m>k jyl 



'-ZjZids}\ 



-^ E E ( 



<7 



^\z,\'ds)'/\ f ^\zi\'dsf'' 



m>k j^l 



^V 



E E / ^i--,r^*+'< E E / 



H^. 



^1^ |2 



zi\ ds 



m>k jyl 

=h31 + -^232 • 






(3.23) 



(3.24) 



(3.25) 



So, we need to estimate the six terms I211, I2i2^ I221, h22, I231, h32- We begin by estimating 
the terms /211 and 1221- We note that 



\U - Uo\la,K Su\ha,K + E (I'^PI ^ll^llo,a;p)l'^P:-^ll,",i<' 



<\ 

r-^ I 



1 



(3.26) 



u\l,a,K + l{a)YjT7\\u\\Q^^^. 



K 



Let u = u — p!— I J^ udx and uq = Rhu. Also note that n — uq = "u — So since Rh preserves 



\^k\ J^k 



constants. It follows from (3.26) and Lemma 3.1 that 



W - uoW^^K = \u- nol? „,^ < \u\l^^^K 



7(a)^^PIIo,c.K <7(a)l^i|?,h,, 



\K 



a,ujK' 



We then have 



hii = I221 = E I"" ~ ""olla,;^ ^ l{oL)Hlxa,u,i ^ l{a)ah{ 



U,U] 



K 
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Next, we estimate the term /232- Estimates for the terms /2i2; -^222 and I231 are similar. 

<!j!ca~K,„fc he 



l232<V Y. e^a^.TT ^ / 



{zil'^ds 



m'>k 



3>l 



< 



We 

max - — 

jiceK„(. he 

mk j^l 

m'>k 



V Yl 



I 1 1 2 II 1 12 



a/c 



Using the trace inequaUty in Lemma 3.1 we have 



W, 



h32<vY Ts¥ ir(^ ^\\^"^\\iK„-.+H\zm\l 



jl'^SK^f. he 

mk j^l 

m'>k 



h,K„ 



dK. 



w, 



+ r/ V max -r^{H ^WzkW^K^, + H\zk\i,h,Kt, 



T7l>fc 



ejiCdK^f. fl. 
mk j'^l 



rri>fc 



max -^{H ^||n-no||o,if^+-H'|n-'Uo|?,h,K„ 



^„ 3lCl>'^mk h, 

rri>fc 



max — - 

mk j-:>l 

m>k 



{H ^\\U - UoWl^Ki, + H\U - Uo\lf„Kk 



aK. 



(3.27) 



(3.28) 



We recall that no is a continuous function on each K G Tff, which implies that 

K - Uo\lxK ^\'^\lxK + \MlxK ^ \u\lhK + l"o|?,i^ 



<Hi,h,K+ Yl (I^pI ^I 

<\u\l,h,K + 1{1)t^A\u\\Iu^^. 



I|2 Mi |2 

u\\o,cVp)\(pp,K\l,K 



(3.29) 



Again let u = n — -rj^-r J^ udx and uq = i?j:/n, and note that u — uq = u — uq {Rh preserves 
constants), we have from (3.29) that 






ILa- ^7(l)l^i|iA<^A'- 



(3.30) 



Using the approximation property of Rh in Lemma 3.2 as well as (3.30), it follows from (3.28) 



that 



We 



l2S2<V E .,--., 7^(^1-1 W.+^7(l)|n|?,,,.,,J 

We 



-^jlCdK^^ he 

"'^mk j>l 

m^k 



+ 77 y^ max -r^(H\u\ 



l^^l^mk he 



Eji^c-^mfc rie 
°^mk j>l 

m'>k 



^,h,i^K^. 



+ Hj{l)\u\U 



^Kt, 



H 



<r/ max max We'y{l)-rah{u,u). 



SKmk '^jl^S'^mk 
m>k jyl 



h 
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Having all the estimates of the terms I21, I22 and I23 together in (3.22), we get 



IT 

/2 < r/ max max We'y{l)—ah{u,u) + -f{a)ah{u,u). 

m>k jyl 



Combing (3.20), the estimates for /i and I2, and Lemma 3.2 we have 



y^ ai{ui,Ui) < (rj max max Wej{l)— + max(j{a),'qP{a)) \ah(u,u) 

I — m>k j^l 

The second assumption of the Schwarz framework, requires a bound for the spectral radius 
p{M) of the N X N matrix M, whose elements Mij are defined in terms of a strengthen Cauchy 
Schwarz inequality: let < Mij < 1 be the minimum values such that 

\ah{ui,Uj)\ < Mijah{ui,Ui)2ahiuj,Uj)2, 
where n, G F^, Uj £ V-", I <i,j < N. 

Note that in our definitions above 

\ah{ui,Uj)\ = 0, if iaieas{dQij) = 0, I < i,j < N, 



where dflij = d^i n dflj. This is because, according to (2.7), all the terms become zero since 
functions Ui and Uj have no common support. For the remaining case we can take Mij = 1. It 
follows at once from Gershgorin's circle theorem that 

p{M) < iVc + 1, 
where Nc is the maximum number of subdomains adjacent to any subdomain. 

The third assumption asks for uj such that 

ah{ui,Ui) < Loai{ui,Ui), Wui G Vj^, , < z < iV. 

Since we use exact bilinear form for the subproblems, cj = 1. 

Our theorem is proved since the analysis of the three assumptions is complete. D 

3.1.2 Overlapping additive Schwarz method 

In this section, we analyze the overlapping version of our additive Schwarz method. As will be 



shown in Theorem 3.2, the condition number bound is not only dependent on the indicators 7(a) 
and /3(a) but also on the partition robustness indicator 7r{a) borrowed from |5j, which describes the 
relationship between the subdomain overlap and the coefficients a. We can control this indicator 
by choosing a suitable overlap. 

Given an overlapping partition {^2'^} with 6i > 0, i = I,- ■ ■, N, let {xi}iLi be the partition of 
unity subordinate to {^^j, c.f., [7], for which the following property holds, 

Xiix) = 1 and Vxi = 0, for all x e n'f = {x e n'i : x ^W- for any j / i}. (3.31) 
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Partition robustness indicator [5]. For a particular partition of unity {xi} subordinate to 
the covering {^[}, let 

1=1 
then the partition robustness indicator is defined as 

7r(a) = inf 7r(a,{x*}), 
{x>}en{{Qa) 

where n({i7^}) denote the set of all the partitions of unity {xi} subordinate to the cover {ri^}. 



Following similarly as the nonoverlapping case, we use Vq as the coarse space with bilinear form 
ao(-, •), the subspaces V^ and the local bilinear forms aj(-, •), 1 < i < N, can be got by taking 



Bi = 0,[ in (3.6) and (3.7) respectively. 



The next lemma gives us an estimate on the boundary layer. 

Lemma 3.3. Let D be a convex domain with diam(D) = 0{H), Td be a family of partitions over 
D, and D = UreT ^ ™^''' diam(r) = 0{h). Let < fj, < H, then for any u £ Vd = IlrcTD^^i'^)^ 
then we have 

I \u\ |o,D^ < C[fiH-'^\ \u\ \Ij^ + fi{n + H)\u\lf,^ij], 

where D^j = {x £ D : dist{x,dD) < /i} denotes the boundary layer of D with width /i, and 



\Ih,d = E.cD l|Vu||2,, + E^.-.cz, /it/ J [u]^ds. 



3>l 



The proof of this lemma can be found in |14j, again the assumption on the convexity of the 
domain D can be dropped, c.f.. 



Theorem 3.2. For all u G V^ , there exists uq G Vq O'^d ui G V^^ , I < i < N, such that 



N 



'^ai{ui,Ui) < CXahiu,u), 



i=0 



Consequently, we have k{T) < A, where k(T) denotes the condition number of the additive Schwarz 



operator T as defined in (3.8), and A is given as 



J-f 

A = 7r(a) max (7(1), /3(1)) max — -^ + max (7(a), ??/3(a)), 

t Oi 

and C is a constant independent of h, H, and a. 

Proof. Again, we use the Schwarz framework to prove this theorem. Like in the nonoverlapping 
case, we estimate the three parameters Cg, p{M) and uj . 

We first estimate Cq . For all u G V^, we may choose mq = Rhu and Ui = Lh{xi{u ~ uq)), where 
Ih is the usual Lagrange interpolation operator. Let w = u — uq and {xi} be the partition of unity 
subordinate to the covering {il'j}, then Ui = //j(xjU)). 
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ai{Ui,Ui) 



For 1 < i < A'^, since Xi = on Fj = di^'j\dQ, we know from (3.7) that 
y / aVui -Vuidx — 2 N^ / {aVui}^[ui]ds 






eC^i 



+ V'^ I ^{Ui)e{Ui)eds. 



eCEi''^ ^^ 



Using the same techniques as in Lemma 3.2, we can write 






3>l 






and 



y^ / Weae(VMi)e • (nj)eneds ^ ^ / aVuj • Viijdx + ^ ^ / -r^{ui)lds 



Consequently, 






——(ui)^ds 






eC=i 



=/l + /2 + Is. 

We begin by estimating term Ii. Since Uj = Ih{xiw), we have 



Ii = 2, / aVui -Vuidx = y^ / a|V/ft,(xjif)| da;. 
Let xj,t = uT J^Xidx, then it fohows from the Bramble-Hilbert Lemma, c.f., [7J, that. 



I _ II <J^-lvxil, yrcn', 



i,5i' 



where Q'^g, = {x £ Q'- : dist[x,dQ.[) < di} denotes the boundary layer of J7^ with width 6i, and 
fl'f = ^'i\^[ g. is the interior part of Q[. 

Since the interpolation operator I^ is stable with respect to the norm 1 1 • | |o oo Q, using the inverse 
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inequality, we have 



a\VIhixi'w)\'^dx 



l«'/'v(4(x, 



w 



IO,T 



\a'/^V{h{x^,rWml, + \\a'/^Vh{{x^ - X^,r)w)\\l 
la^^^'^wWl, + Wa^/^hiiXt - Xi,r)u^)||o,oo,r 

\a^^'^Vw\\l + K^aWxi - Xi,r||o,oo,rlkllL, 



which imphes that 



/ 



a.|V4(,.„)|^d. < ||oV^V„||S,, + JHVv.nkllt. ^^^^A' 



Adding this estimate across all the fine elements r C 0^, we have 

h=y~] aVIhixiw) ■VIh{Xiw)dx <\\a\Vxi\'^\\o,oo,n\\w\\on' +l'^liaO'- 



rcn', 



By the definition of vr(a), and using the estimate in Lemma 3.3 with fi = 6i and D = il.'-, we 
have 

1 






Sf 



Hi 



^ I \Hi sr^ I |2 , / \^i\ \2 I I |2 



(3.32) 



where in the last inequality we have used the approximation property of Rh from Lemma 3.2 
Next, we estimate the term l2- We first note that 



l\\2 



INIo,e,, = H- n^|^,e,, = Ihixiw^) - h{XiW% 

= IhiXiiw^ -W^))\o,e,i < \Wj -Wl\l,e,i = l[Hlo,e,,> 



(3.33) 



where in the above equality we have used |xj| — ^i 1 — ^ — -^) ^^^ ^^^ stability property of Ih 
w.r.t. the norm || • ||o,oo,c- 



It follows from (3.33) that 



w. 






[uifda-^-q Y 



W, 



e r„„i2 



w] ds. 



J>1 
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Finally, for the term I3, since, by definition, Ui = Ih{xi'w) and Xi = 1 on Hj, we have 



'^ = ^T.l^(-^)"ds<,Y.l 



We 2 , 

—Wgds. 



Combining the estimates of /i , I2 and I3, 1 < i < N, we have 

N N N 

3>l 



Uj] ds 



+ rj 

i=l eC; 

^ H ^ H ^ 

<^7r(a)^ E l^l?A-K+E^(«)7^l^l?Af^:+El^lLo: (3-34) 

— w^as 



Hi ^ Hi 

<7r(a)max— ^a/,(M,u) + ^^71(0)— ^|u - iJ^^li ^n' +ah{w,w)- 



Note that, 

^ Hi ^ Hi ^ Hi 

^T^{a)^\u- Rhu\\^^, <^7r(a)-^|n|2 , +^7r(Q)-^|i2Hu|2 ,. (3.35) 

^ Oj ' ' « ^ dj • • » ^ Oi ' ' » 

1 = 1 2=1 J = l 

It follows from Lemma 13.21 that 

N 



^ \RHu\\h,nr ^ max (7(1), /3(l))a/,(u, u) 

i=l 



which, together with (3.35), implies that 

X^ vr(a)-T^|M - Rhu\Ij^q> < 7r(a) max — ^ max (7(1), (3{l))ah{u, u). 

■ -, Oi ■ ■ I I Oi 

2 = 1 

Since uq = Rhu, using the stability property of Rhu in Lemma [3^ we have 

a/i(it;, w) =ah{u -uo,u- uq) < ah{u, u) + ah{uo,uo) 
< max (7(0), ?7/3(a))a/i(n, u). 



Thus, by (3.34), we have 



H 



'^ai{ui,Ui) < 7r(a)max— ^max(7(l),/3(l)) + max (7(a), ?7/3(a))aft(u,n). 



j=i 
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Note that ao{RHU, Rhu) < max ['y{a),r]P{a))ah{u,u), which imphes that 
Cq ^ 7r(")max(7(l),/3(l)) max — ^ + max (7(a), r/ /3(a)). 

i Oi 

The other two parameters p{M) and w are estimated in the same way as before. D 

3.2 Oscillatory boundary conditions 

In this section, we foUow the notations used in [5], and give an exphcit bound for the indicator 7(a) 
with a shghtly different proof. We show that if the high-conductivity region crosses the boundaries 
of coarse grid blocks, the coarse basis function with linear boundary condition fails to give a robust 
bound for the condition number. The coarse basis functions with oscillatory boundary condition, 
on the other hand, yield a robust method. First, for each K € 7}/, let p > 1 be an arbitrary 
constant, define the set 

K{p) := {x £ K, a{x) > p}. 

Since a{x) is piecewise constant with respect to Th, K{p) is a union of fine grid elements. Let the 
region 

K{p) = K\p)UK''{p) 

be associated with each K G Th, where the set K^{p) contains the components of K{p) whose 
closure touches dK and K (p) contains all the interior components of K{p). The term e{p,K), 
representing the distance between K^{p) and K^{p), be defined as 

e{p,K) = dist{K'{p),K^{p)). 

Assumption 3.1. (1) K (p) and K {p) should be well-separated, i.e., 

e{p,K)>h. 

(2) K (p) can be written as a union K (p) = '^i=iKi{p), where the components Kf{p) are 
simply connected and pairwise disjoint, and a(x) is a constant on the closure of each component, 

i-e., 

a{x) = ai, for all x G Kj^{p). 

(3) Let Tf{p) = KP{p) n dK be the boundary part of dK which locates in KP{p), and that 
T^{p) = "^i^i Tf{p). For all K G Th, let E be any edge of dK, we require that 

\E\T''{p)\>Hk, 

which means that the high conductivity field does not cover too much of dK, \/K G Th- 

Note that, our assumptions are weaker than those [5|, we do not need the coefficient q(x) to 
be continuous across the coarse grid boundaries. Next we give the explicit bound for the indicator 



7(a) in Theorem 3.3 with a different proof than in [5j. 



Theorem 3.3. Let Assumption 3.1 hold true for each K G Th, then an upper hound for the 
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indicator ^ (a), with ipp^oK being the oscillatory boundary condition, is given by 



7(a) < P 



H 



Proof. The key idea is to partition K{p) into two parts Ki{p) and K2{p), c.f., Figure 1, build a 
special function 6 € S^{K) whose bound can be estimated. 




K2{P) 



Ki(p) 



Figure 1: The areas surrounded by dashed lines are K^{p) and K^{p) respectively. Ki(p) and 
K2{p) are separated by dotted lines. 

Let K S Th, and tppOK be the oscillatory boundary condition, (pp k is the multiscale basis 



function which is built through the discrete harmonic extension, c.f., (3.3). Accordingly, 



t>p,K\l,a 



K < |^|i,a,i<', for all ^ G 5 {K) which satisfies 6\qk = 4>p,K\dK- 



We only need to construct a function 9 for which we can estimate its | • \i,a,K norm. We define 
the function 6 explicitly by its values at the nodes oi J\f^{K), where M^{K) denotes the set of fine 
mesh vertexes in K. Obviously, the function 6 is contained in S (K). We begin by constructing 
e on each K^{p). By Assumption 3.1, K^{p) = Ef=i^f (p) and that r^{p) = Y.Li'^fip) ^it^ 
r^ = KP{p)ndK being the part of dK locating in Kj^{p). Next, we define the values of 6 on Tf{p), 
then extend to Kj^{p), \/ 1 < I < L. For simplicity, we assume that F^ lies only in the interior of 
dK, i.e., rf does not touch any vertex of K. For a fixed /, we may choose local coordinate system 
(xi,X2) and some hi, 62 > such that K C {(xi,X2) : X2 > 0}, F^ = {(xi,0) : xi G [bi, 62]}- Define 
9{xi,X2) = TpixijO) for all (xi,X2) G Kj^{p). After defining all the values of 6 on K^{p), its values 
at the nodes of K\K^{p) for which we set 0. 
Define 

Klip) = {xeK : d\st{x,dK) < d}, 

where d = max „g, , dist{x,dK). Similarly, i^2(p) = K\Ki{p) denote the area of K which is not 

contained in Ki{p). Clearly, K^{p) C Ki{p) and K\p) C i^2(p) hold by their definitions. 
We give an upper bound for \9W ^ j^. We note. 



l^ll,a,i^ < l^ll,a,Xi 



(p) 



\l,»,K2{p) 



h + h 
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For the term Ji, we have 

|2 |/)|2 



h - \0\l,a,Ki{p) - \^\l,a,KB{p) + \^\l,a,Ki(p)\KB{p) " -^11 +-^12- 

From (3.2) and the nodal values defined above we can see that (g^)^+(g^)^ = o.""^ / { Jy{ci^ )^^ ds) 
on each r C K^{p), where T n K^{p) 7^ is the coarse edge of K. Hence 

In= Y. j <ij^f + ^^)fdxidx2= Y. [c^-'/i [(.c.'^r'ds)^dxidx2. 

TeKB(p)'^'^ ^ y rdKB^p)-''^ ■''^ 

Let Ti = 1r\K^{p) and T2 = T\K^{p), since (q^)^^ > p^^ on T2, it follows from Assumption 
3.1 (3) that 

! {a')-^ds= [ {a'^)-'ds+ [ (a'^y^dsy [ {a^y^ds > p-'\T2\ > p'^Hk, 
Jr JTi JT2 Jt2 

which implies that 

Iu= Y [c^-'/{[ic^^r'dsYdxidx2< Yl P''h\/P''Hl<p. 



For the term lu, let K^ip) = Ki{p)\K^{p), and 

K^iph) = {rC Ks{p) : t n dK^ip) / 0} 

be the boundary layer of K^{p) with width h. Note that, by definition, we have V0 = on 
K3{p)\K3{ph) = {t C K^{p) : fr\dK-i{p) = 0}, and since a < p on K^ip), we get 

TGi^3{p) r(iKz(ph) reK^iph) 

where in the last inequality we have used the fact that \0\ij- ^ 1 since all the nodal values of 9 
defined on r are between and 1. 
Finally, for the term I2, define 

K^ip) = {tC K2{p) ■.TndK2{p) / 0}, 

which is the boundary layer of width h of K2{p), and 

Kiip) = K2{p)\K^ip), 

which is the interior part of K2 [p) ■ 

By the definition of 9, we have 9{x) = on i^|(p). Since a < p on i^^(/)), we get 

h = \9\l^,K,^,) = E \^\Ur<P E \^\lr^Pj- 
rCKB(p) rCKB(p) 

The theorem is proved by combining the upper bounds for /i and I2 together. D 
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3.3 Upper bound of DG indicator 

In this section, we estimate the term /3(a). Note that, if the coefficient field is continuous across 
the coarse grid boundaries, we have /3^(a) = 0, thus in this case, what we need to focus on is the 
estimate of the term /3 (a). However, in reality, the coefficient field can be discontinuous across 
the coarse grid boundaries. In both cases, suppose the coefficient of the conductivity field satisfies 
We ^ p, where p > 1 is a given number, on the edges of fine mesh triangles which intersect with 
the coarse grid boundaries. Then the indicator /3(a) yield a bound which is independent of the 
high contrast in the coefficients. 

Theorem 3.4. For all K G Tff, e C dK, let the coefficient field he discontinuous across the coarse 
grid boundaries. If there exists p > such that We ^ p, then 

/3(a) <pf. 



Proof. By definition, c.f., (3.9), we have 



B (a) = max max > ( / ^-^ I [(/>ol Pds) 

< -,„ -. E (/fiwi^.,) 

a^mfeCn ^p&JsfHi^K^) ^ J he 

m>k jTjy <=]lCaKmk 

^flii- , Je "-e 



H 

3>l 



~ .™ax^ max 7 ^ (^ / —i-as) ;;; p 



xpeaK^f. J 



using the fact that \[(t>p\\ ^ 1- The same result holds for the estimate on /3^{a). Hence, since 
/3(a) = max (/3^(a),/3^(a)), the theorem is proved. D 

4 Numerical Experiments 



In this section, we present our numerical results, where we solve the equation (2.1, 2.2) with g = 
on the square domain ri = [0, 1]^. We run the preconditioned conjugate gradient method until the 
I2 norm of the residual is reduced by a factor of 10^. 

In the numerical experiments, subdomains {Qi}^^ are all square shaped, and each subdomain 
consists of two coarse triangles. In each of our numerical experiments below, we consider the 
performance of our two methods and compare with the method in [5J. The results of our method 
are shown in the columns under the heading "Discontinuous Galerkin", while the results from [5j 
are shown under the heading " Continuous Galerkin" . We also use different basis functions for the 
coarse space, 4>p, (j)p ' , and (j)p ' , that is, the piecewise linear basis function, the multiscale 
basis function with linear boundary condition, and the multiscale basis function with oscillatory 
boundary condition, respectively. We choose the same penalty parameter ry for both the fine and 
coarse bilinear forms. 

Example 4.1. We begin with our first example here, see Figure 2, where H > 8h and a(x) denotes 
a 'binary' medium with a(x) = a on a square area inclusion lying in the middle of each coarse grid 
element K G Th, and at a distance of H/8 both from the horizontal and the vertical edge of K, 
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Figure 2: The coefficient a{x) corresponding to the binary domain of Example 4.1. 

and a(x) = 1 in the rest of the domain. We study the behavior of the preconditioners as a — t- oo. 
For the discontinuous Galerkin formulation, we first consider the nonoverlapping method. We 
note that maxx^j. maXec^K^j. We < 1, and that the DG indicator /3(a) = because the coefficient 

JTi>fc 

is continuous across the c oars e grid boundaries as well as the zero Dirichlet boundary condition. 

implies that 7^^' (a) < f . 



Taking p = 1 in Theorem 



3.3 



3.1 



that 



It then follows from Theorem 
the nonoverlapping method has a bound which is independent of the jumps. For the overlapping 
method, we have 7r(a) = (5?||a|VxiP||L°°(Q) ^ li thus in this case the multiscale basis function with 
linear boundary condition will yield a robust bound. We note that, in this example, the multiscale 
basis function with oscillatory boundary condition is the same as the one with linear boundary 
condition. 

The numerical results in Table 1 show that, both the nonoverlapping and overlapping method 
with multiscale coarse basis functions with linear boundary conditions are robust as predicted by 
the theory. The overlapping method [5j with the same multiscale coarse basis functions produces 
almost the same condition number estimates, however, for the linear coarsening, the results in 
Table 1 show a loss of robustness of the overlapping Schwarz method as a goes from lO" to 10^. 



Table 1: Condition number estimates of the Schwarz methods on Example 4.1 with h = 1/128, 
H = 8h, 6 = 2/1, and 7/ = 4. 







Discontinuous Galerkin 


Continuous 


Galerkin 




Nonoverlapping 


Overlapping 




Overlapping 


a 


MS,L 




MS,L 


MS,Lr N 


MS,L 


L 


10° 


26.79 




6.43 


1.00 


5.64 


5.6 


102 


21.61 




6.44 


1.42 


5.79 


58.6 


10^ 


21.49 




6.81 


1.44 


5.80 


358.3 


10^ 


21.49 




6.81 


1.44 


5.80 


378.7 



In our next experiments, we study the behavior with different penalty terms. As we can see 
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from Table 2, the condition number estimate for the nonoverlapping method grows linearly with 
the penalty parameter rj. However, it is almost constant for the overlapping method, which suggests 
that the condition number bound do not depend on the penalty parameter. The results in Table 2 



is in agreement with Theorem 3.1 and Theorem 3.2 



Table 2: Discontinuous Galerkin formulation on Example 4.1. Condition number estimates of the 
Schwarz method with h = 1/128, H = 8h and 5 = 2/i(only for overlapping). 





Nonoverlapping 


method 


Overlapping 


method 


a 


7? = 5 


r] = 10 


7]= 100 


r? = 5 


r/= 10 


7/ = 100 


10° 


33.14 


64.88 


635.2 


6.47 


6.56 


6.88 


10^ 


26.55 


51.09 


487.3 


6.73 


6.87 


6.96 


10^ 


26.39 


50.74 


483.4 


6.83 


6.89 


6.99 


10^ 


26.40 


50.73 


483.4 


6.83 


6.89 


6.99 



Remark 4.1. Note that, in our numerical experiments, we test the model problem with zero 
Dirichlet boundary condition, thus all the degrees of freedom are inside domain 0, which means 



that our DG indicator /3(a) defined in (3.11) is only for dKmk C $7 but not for the coarse grid 
boundaries on dVl. Moreover, in Example 4.1, since the coefficient of the conductivity field is 
continuous across coarse grid boundaries, we have /3(q) = due to the fact that the jump of the 



coarse basis functions corresponding to the same coarse node is equal to 0. From Theorem 3.1 we 
know that the condition number bound for the nonoverlapping case will be 

r/max max We7(l)— + 7(a), 

m>k jyl 



taking p = 1 in Theorem 3.3 we have 7(a) ^ -^, thus theoretically speaking, the condition number 



bound will grow linearly with the penalty parameter 7]. While for the overlapping case. From 



Theorem 3.2, we know that the condition number bound for the overlapping case will be 



J-f 

7r(a) max(7(l), /3(1)) max — — h 7(a) 

i di 



taking p = 1 in Theorem 



3.3 



H 



implies 7(a) < ^, thus in this case, our condition number bound 
will keeps unchanged independent of the penalty parameter r/. Both of the two cases are matched 
exactly by our numerical experiments in Table 2. 



Example 4.2. In this example, the high conductivity field is no longer contained inside, it touches 
the coarse grid boundaries, see Figure 3 and Figure 4. In Figure 3, The conductivity field a{x) 
corresponds to a binary medium with background a{x) being equal to one, and a = S on the 
channels with an area equal to 2h x 4/i each. In Figure 4, the high conductivity channels with 
diameter 2/i x 2h are located only on one side of the coarse grid boundaries. 

We first compare the behavior of our nonoverlapping method on the two cases shown in Figure 3 
and Figure 4. For the the conductivity fi eld g iven in Figure 3, 7(a) < -^ since max^ laaXecdK We 
c^i li'^) ^ 'h (taking p = 1 in Theorem 
condition number bound will grow as a 



3.3), and /3(a) = 0. Thus Theorem 3.1 predicts that the 



00. However, for the conductivity field in Figure 4, we 
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Figure 3: Example 4.2 with a{x) crossing 
the coarse grid boundaries, a{x) = a on the 
channels and a{x) = 1 otherwise. 



Figure 4: Example 4.2 with a{x) touching 
the coarse grid boundaries, a{x) = S on the 
channels and a{x) = 1 otherwise. 



have mayiK niaXecdK We < l- 



P 



7(a) 



1 in Theorem 3.4). 



H 

h 



(taking /> = 1 in Theorem 3.3), and f3{a 



< 



f (taking 

Theorem 3.1 predicts that the bound is robust in this case. For the 
overlapping method, it is easy to see that 7r(a) < maxj (5?||aV|xiP||Loo(^) ~ ^«V^^ ~ 1- It follows 
from Theorem |3.2| that the condition number bound will be independent of the jumps for both 
the conductivity field in Figure 3 and in Figure 4. The numerical results in Table 3 confirm these 
results. 



Table 3: Discontinuous Galerkin formulation on Example 4.2. Condition number estimates for the 
Schwarz method with h = 1/128, H = 8h, 6 = Ah, and the penalty parameter r] = A. 



Nonoverlapping Overlapping 
a Fig 3 Fig 4 Fig 3 Fig 4 



10° 2.68e+l 26.79 6.23 6.24 

5.00e+2 28.91 7.17 7.15 

4.76e+4 28.95 7.23 7.17 

4.76e+6 28.96 7.23 7.17 



10^ 
10^ 
10^ 



Example 4.3. [Gaussian random field] In this example, we test our method on a more realistic 
model. The coefficient a is a realisation of a log-normal random field (see Figure 5), i.e., logQ(x) is 
a realisation of a homogeneous isotropic Gaussian random field with spherical covariance function, 
mean 0. This is a commonly studied model in the multiscale area in many literatures [5l|6]. There 
are some evidences from field data that this gives a reasonable presentation of reality in certain 
cases |15| ll6j. There are many good ways to generate such random fields, we simply follow the way 
in [19] which used FFT (Fast Fourier Transformation) method \\.1\ [T8] . The spherical covariance 
function has a parameter 0, a bigger 9 increases the correlation of the random field. 
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Figure 5: 256 x 256 Gaussian random field with covariance parameter 6 = 50. 

Table 4: Condition numbers (iteration numbers) of Schwarz methods on Example 4.3 with h 
1/128, H = 8h,6 = 2h, and rj = A. 





Discontinuous Galerkin 


Continuous Galerkin 




Nonoverlapping 


Overlapping 


Overlapping 


High Contrast 


MS,OSC 


MS,OSC 


MS,OSC 


8.55e+3 


30.72 (52) 


7.20 (16) 


6.65 (15) 


1.75e+5 


32.68 (57) 


7.85 (18) 


7.12 (16) 


7.32e+7 


36.40 (67) 


9.86 (21) 


8.45 (19) 


1.49e+9 


37.90 (73) 


11.27 (22) 


9.38 (20) 


6.02e+ll 


45.94 (86) 


14.77 (25) 


11.73 (22) 


1.28e+13 


51.61 (96) 


16.80 (26) 


13.28 (25) 



In the numerical experiments, we compared the condition number estimates (iteration numbers) 
of both the nonoverlapping and overlapping domain decomposition method with the method [5j. 



5 Conclusions 

In this paper, we present two two level additive Schwarz domain decomposition methods for the 
multiscale second order elliptic problem. The work here is an extension of related works, c.f., [H |5], 
to a discontinuous Galerkin formulation, c.f., [13]. For the nonoverlapping case, we show that if the 
conductivity field do not cross the subdomain boundaries, our nonoverlapping method is robust. 
However, in some cases, when the conductivity field cross the subdomain boundaries, we show 



in Theorem 3.1 that our condition number estimate will be dependent on We, i.e., the weighted 
average of the coefficient over the edge e. It is now the topic of further investigation to see if we can 



get rid of this dependence on We- For the overlapping method, as seen in Theorem 3.2 we show 
that our method is robust under Assumption 3.1. We allow the discontinuity of the coefficients 
across the coarse grid boundaries. 
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